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ABSTRACT 


This report consists of a series of graphs whose purpose is to 
indicate the azimuth and elevation of the sun for any day in the year, 
for anytime of day, at Wallops Island Range, the Churchill Research 
Range, or the White Sands Missile Range. An indication is also 
made of the azimuth and elevat.on of the geomagnetic vector at each 
of the launch sites. The report includes instructions for use of the 
graphs , and a section explaining the origin and the computation of 
the coordinates for the graphs. 
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POSITION OF SUN AND MAGNETIC VECTOR FOR THREE LAUNCH SITES 
(Wallops Island, White Sands, and Fort Churchill) 


INTRODUCTION 

One of the documentation requirements of many rocket-borne experiments 
is that of recording the vehicle attitude at the time an observation is made. Data 
from the solar sensor and the magnetometer, probably the two most widely used 
aspect sensors in the sounding rocket field, must be used together with the 
knowledge of the position of the sun and magnetic vector in determining the atti- 
tude of sounding rockets. 

The purpose of this document is to provide sufficient information so that a 
user of sounding rocketry will be able to launch the vehicle at a desirable time 
of day, as the time affects the direction of the sun in relation to the direction 
of launch. 


FORMAT OF SOLAR DIRECTION DATA 


Solar direction data are plotted in graphs, each containing a family of curves 
showing the direction or location of the sun as a function of the launch site, as 
well as the day of the year, and the time of day. 


HOW THE GRAPHS ARE CONSTRUCTED 

The curves are plotted according to a coordinate system in which the ab- 
scissas are in azimuth degrees, and the ordinates are in degrees of elevation 
with reference to the horizon. Each curve represents the apparent path taken 
by the sun on one particular day of the year. There are 18 curves for one year, 
with 20-day intervals between the days represented by adjacent curves. Two 
graphs are made for each site: one for January through June, the other for July 
through December. 

Curves for each of the selected days are generated by plotting points at one- 
hour intervals over a 24-hour period on a scale that places 180 degrees azimuth, 
or south, at the center of the horizontal scale, and zero degrees of elevation, or 
horizontal, at the center of the vertical scale. Then, all of the hourly points 
representing a given day are connected by a curve. After all.of the curves for a 
given graph have been drawn in a similar manner, all of the points representing 
each hour are connected by a curve. Each of the day curves is labeled accord- 
ing to the day for which that curve is valid. Each curve representing the hour 
is also identified. 
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On each graph, another point is plotted to show the azimuth and elevation of 
ihe magnetic vector at the rocket range for which the graph is valid. The com- 
pleted graphs are shown in Appendix Figures A-l through A-6. 


HOW TO USE THE GRAPHS 


To find the direction of the sun on one of the days for which a curve is plot- 
ted, at an integral hour point, find the coordinates of the encircled point. If the 
desired time includes a fraction of an hour, move a proportional distance along 
the line connecting two hour-points and read the coordinates of the point thus 
located, on the horizontal and on the vertical scales of the graph. 

To locate the sun at an integral hour point , but on a day other than one of the 
days represented by a curve, move a proportional distance between two adjacent 
day -curves along the line connecting each day-curve for the desired hour. When 
that point on the line is located, find the angles for the point on the horizontal 
and vertical scales. Figure 1 illustrates interpolation as explained here. 

The azimuth and elevation of the magnetic field are shown on each of the 
graphs. Table I gives the directions of the magnetic fields of the three launch 
sites, at different heights. 


DETERMINING THE LAUNCH WINDOW 

Two conditions should be observed in choosing the time and the direction of 
launch which are compatible with the use of solar and magnetic sensors. 


1. The angle between the direction of the sun and the longitudinal, or Z, 
axis of the rocket, should be great enough that the sun will be within 
the field of view of the side-looking solar sensor. 

2. The angle between the direction of the sun, and that of the magnetic 
vector, should also be of sufficient magnitude to produce a well-defined 
intersection in the pattern indicating the angles of the sun sensor and of 
the magnetic sensor. This will allow an accurate determination of the 
attitude of the rocket. 
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DEGREES OF ELEVATION 



POINT 






AZIMUTH 

ELEVATION 

A 

= SUN 

LOCATION AT 

1000 

ON APRIL 

17 

144.3 

59.9 

B 

= SUN LOCATION AT 

1033 

ON APRIL 

17 

154.3 

66.2 

C 

= SUN 

LOCATION AT 

1000 

ON APRIL 

5 

149.4 

53.2 

D 

= SUN 

LOCATION AT 

1033 

ON APRIL 

5 

158 

58.7 


Figure 1. Interpolation on Graphs 
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TABLE I 

MAGNETIC ELEVATION ANGLES (1970) 


WALLOPS BLAND, VIRGINIA 
37.8’ NORTH LATITUDE, 75.5* WEST LONGITUDE 

ALTITUDE 

ELEVATION 

AZIMUTH 

1 

TOTAL FIELD 

KILOMETERS 

DEGREES 

DEGREES 

MILLIGAUSS 

0 

69.7 

171.9 

556.4 

100 

69.6 

172.1 

528.4 

200 

69.5 

172.3 

502.4 

300 

69.4 

172.5 

*78. 0 

FORT CHURCHILL, CANADA 
58.75° NORTH LATITUDE, 94.0° WEST LONGITUDE 

ALTITUDE 

ELEVATION 

AZIMUTH 

TOTAL FIELD 

KILOMETERS 

DEGREES 

DEGREES 

MILLIGAUSS 

0 

83.4 

184.1 

615.0 

100 

83.2 

184.4 

585.2 

200 

83.0 

184.5 

557.3 

300 

J 

82.8 

184.7 

531.2 

WHITE SANDS, 
32.38° NORTH LATITUDE, 

NEW MEXICO 
106.5° WEST LONGITUDE 

ALTITUDE 

ELEVATION 

AZIMUTH 

TOTAL FIELD 

KILOMETERS 

DEGREES 

DEGREES 

MILLIGAUSS 

0 

60.5 


520.7 

100 

60.6 


494.9 

200 

60.6 

191.6 

470.7 

300 

60.6 

191.5 

448.2 


NOTE: The South-Seeking Direction, or Tail of the vector, is shown, since it 
has a positive elevation angle. 
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The time of day of the launch is under the control of the experimenters, un- 
less the time is dictated by the requirements of the experiment, or by range 
scheduling problems. Any choice of launch azimuth angles is severely limited 
by range safety requirements, but the launch azimuth may usually be varied over 
a email angle. Consequently, it is usually possible to sufficiently change either 
the launch time, or the launch angle, so that the relationship between the direc- 
tion of launch and the direction of the sun will ensure adequate performance of the 
solar and magnetic sensors. 


DERIVATION OF GRAPH COORDINATES 


THEORY 

The computations involved in compiling the tables of coordinates are, es- 
entially those of the science of navigation. This science enables a traveler to 
determine his location on the surface of the earth by observation of the position 
of certain celestial bodies, at a precise time. There are two publications which 
are usually used to interpret the data obtained by such observations: 

1. The Air Almanac (Reference 1), gives the directions of the sun, and a 
number of other celestial bodies, from the center of the earth, as a 
function of time, and expressed in terms of a set of universal 
coordinates. 

2. American Practical Navigator (Reference 2) is a reference work on 
spherical trigonometry as it relates to terrestrial navigation. Use of 
the methods contained in Bowditch translates the direction of the sun, 
from the coordinates of the Air Almanac, to coordinates expressed with 
respect to a specific location on the earth's surface. 


DETERMINING THE COORDINATES 

Coordinates for the graphs of Figures 1 through 6 were obtained by solving 
two spherical trigonometrical equations taken from Bowditch pages 28 and 31. 
These two equations were chosen because they were considered especially suited 
to computer methods. For the angle of elevation, 

Sin E = Sin L Sin D + Cos L Cos D Cos T , ( 1 ] 
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and for the azimuth angle, 

Sin (D-L)Cos 2 T + Sin (D+L) Sin 2 T 

Cos A * 1 - 1 2 ) 

Cos E 

E ® Angle of elevation (-90° to +90*), degrees above and below the horizon. 

A = Angle of azimuth (0 # to 360°), measured clockwise from north. 

D a Declination of the sun, 90 degrees north(+) to 90 degrees south(-) of 
the equator. 

L = Latitude of the site, 90 degrees north(+) to 90 degrees south(-), refer- 
enced to the equator. 

T = Ix>cal hour angle of the sun, 180 degrees east(+) to 180 degrees west(-), 
referenced to the launch site. 

Latitude and longitude of each of the three sites are contained in Table I. 

Declination and Greenwich hour angle of the sun are obtained from the Air 

Almanac. 

The local hour angle is found by algebraically adding the Greenwich hour 
angle to the longitude of the site. 

COMPUTER PROC 1AM. Figure 2 illustrates a computer program to take 
the two site coordinates and the two available solar coordinates as inputs, and 
give the desired outputs. It is written in A Programming Language (APL), the 
language of the APL/360 computer system (Reference 3). 

The program defines a function which, when provided with numerical values 
for each of the parameters: latitude and longitude of the site, and declination 
and Greenwich hour angle of the sun, will give numerical values for the angles 
of elevation and azimuth of the sun. 

EXPLANATION OF THE PROGRAM. Steps 1 through 8 of the program 
illustrated in Figure 2 are for the assigning of numerical values to the four 
variables. 
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Steps 9 through 11 convert degrees to radians. 

Steps 12 through 14 compute the local hour angle in radians. 

Steps 15 through 17 solve equation 1, and convert the answer to degrees, 
with the label of ALT. 

Steps 18 through 22 solve equation 2. Its answer : in degrees, is labeled 
AZT. 

Steps 23 and 24 give the answers, in elevation and azimuth angles. 

USE OF PROGRAM. The program can be stored, using the key word, SUN. 
Figure 3 contains an example of how the program can be used to solve the two 
equations for each pair of coordinates for the graphs. 


ERRORS CAUSED BY LEAP YEAR CYCLE 

Because the length of a year is not an integral multiple of the length of a day, 
the leap year system was developed to keep the long-term time relationship be- 
tween the day and the year nearly constant. However, over a four-year period, 
there is a slight variation of the time relationship between any given day and the 
solar phenomena, the elevation and azimuth angles, but the greatest variation, 
that of the elevation is not considered great enough to be significant. Table II 
shows some typical five-year comparisons for a specific day and hour. The 
maximum difference in elevation over the period 1962 through 1966, is found to 
be 0.29 degrees, or 17.4 minutes. 
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VSUNll}} V 
7 SUN 

ClJ » ENTER DECLINATION *00 TO 90 • 

[ 2 J DEC * □ 

[ 3 J • ENTER GREENWICH HOUR ANGLE 0 TO 360 ' 

[4j GHA+ n 

C 5 J ' ENTER LONGITUDE "180 TO 180 • 

C 0 J L0N + U 

C 7 J • ENTER LATITUDE "90 TO 90 * 

CO] LAT+ft 

[9] RA D+ 57,2957795 
[10J D+DEC*RAD 
[llj L+LAT*RAD 
[12] TTT+GHA+LON 

[ 13 J LHA + i TT?*(TTTi 100) ) ♦ ( “ 3 &Q+TTT ) * ( TTT> 1 0 0 ) 

[14] T+LHA *RAD 

[ 1 6 J //«-( ( 20/) ) x ( 2oL ) x ( ?oT) )♦( 10D)*10L 

[10] /l-’lotf 
[ 17 ] ALT+A*RAD 

[10J ZG*-( << lo(D-L)) x((2o(7’* 2 ))*?))♦( ( lo( />♦/,) )x((lo(7« 2) )*?.)))* 20 / 

[19] z«*"2o;$s 

[ 20 ] At1+-RAD*Z 

[21] AZH+-3C>Q-AH 

[ 22 J /lZr«-( ( ( T < 0 )* Atl ) + ( ( r^0 ) x/l 2// ) ) 

[23] ' ELEVATION ' \ALT 

[ 24 ] • AZIMUTH ' lAZT 
V 


Figure 2. APC/360 Program for Computing Elevation and Azimuth of Sun 


SUN 

ENTER DECLINATION *90 TO 90 

Us 

* 21 . 8 

ENTER GREENWICH HOUR ANGLE 0 TO 360 

[]: 

178 

ENTER LONGITUDE "180 TO 180 

□ j 

*94 

ENTER LATITUDE *90 TO 90 

□ : 

58.8 

ELEVATION * 15.51 
AZIMUTH 253.4 

Figure 3. Example of Use of Program 
Shown in Figure 2 
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TABLE n 

COMPARISONS OF SOLAR COORDINATES, SHOWING 
VARIATIONS DURING A LEAP YEAR CYCLE* 


DAY 

YEAR 

GREENWICH 
HOUR ANGLE 
(DEGREES) 

LOCAL HOUR 
ANGLE 
DEGREES 

ELEVATION 

DEGREES 

AZIMUTH 

DEGREES 

Jan. 11 

1962 

88.00 

-6.00 

9.247 

174.4 

Jan. 11 

■ ■■ ■ — — 

1963 

88.03 

-5.97 

9.219 

174.4 

Jan. 11 

1964 

88.07 

-5.93 

9.181 

174.4 

Jan. 11 

1965 

87.97 

-6.03 

9.295 

174.3 

Jan. 11 

1966 

88.00 

-6.00 

9.267 

174.4 

April 1 

1962 

89.01 

-4.99 

35.63 

173.9 

April 1 

1963 

89.00 

-5.00 

35.54 

173.9 

April 1 

1964 

89.05 

-4.95 

35.83 

173.9 

April 1 

1965 

89.03 

-4.97 

35.73 

173.9 

April 1 

1966 

89.01 

-4.99 

35.63 

173.9 

July 10 

1962 

88.70 

-5.30 

53.23 

171.8 

July 10 

1963 

88.70 

-5.30 

53.38 

171.8 

July 10 

1964 

88.67 

-5.33 

53.17 

171.7 

July 10 

1965 

88.68 

-5.32 

53.20 

171.8 

July 10 

1966 

88.68 

-5.32 

53.23 

171.8 

Oct. 18 

1962 

93.70 

-0.30 

21.57 

179.7 

Oct. 18 

1963 

93.68 

-0.32 

21.65 

179. 1 

Oct. 18 

1964 

93.72 

-0.28 

21.37 

179. Y 

Oct. 18 

1965 

93.70 

-0.30 

21.47 

179.7 

Oct. 18 

1966 

93.70 

-0.30 

21.55 

179.7 


•The time of day for each entry is 1800, or 6:00 P.M. local time. 
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Figure A-l. Positions of Sun, January-June, Wallops Island Station 
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Figure A-4. Positions of Sun, July* December, Fort ^hurchill Station 
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Figure A-5. Positions of Sun, Januory-June, White Sands Station 
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Figure A-6. Positions of Sun, July-January, White Sands Station 
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